Coherent Destruction of Photon Emission from a Single Molecule Source 
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The behavior of a single molecule driven simultaneously by a laser and by an electric radio 
frequency field is investigated using a non-Hermitian Hamiltonian approach. Employing the renor- 
malization group method for differential equations we calculate the average waiting time for the first 
photon emission event to occur, and determine the conditions for the suppression and enhancement 
of photon emission. An abrupt transition from localization-like behavior to derealization behavior 
is found. 
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Periodic driving of a two level system can significantly 
alter its dynamics. Two renowned examples are: dynam- 
ical localization of ultra-cold atoms in magneto-optical 
traps ]], and coherent destruction of tunneling (CDT) 
for a particle (e.g. an electron) in a bistable system Q. 
In the latter example, a sinusoidal external driving field 
localizes the electron in one of the wells of a double well 
potential. This is done via particular choice of parame- 
ters of the driving field Q. Such a coherent localization 
behavior was recently addressed in Ref. 0] in the context 
of single molecule spectroscopy 0. 

When a molecule interacts with a continuous wave 
laser field and a much slower radio frequency (rf) field, 
one can control the emission of the molecule by the ap- 
propriate selection of the rf field parameters (see j.3]). 
Coherent localization in this case manifests itself as the 
divergence of average waiting time for a single photon 
emission, even when the laser frequency is in resonance 
with the electronic absorption frequency of the molecule. 
The rf field effectively localizes the wave function of the 
molecule in its ground state (i.e. for certain special pa- 
rameters of the rf driving field) . Note that control of the 
emission by single molecules using rf fields has gained re- 
cent interest since it can be used for single photon control 
and hence, possibly, for quantum cryptology and quan- 
tum computing 0.01- 

In this Letter we present the criterion for the destruc- 
tion and enhancement of emission from a single molecule 
source interacting with an on-resonance continuous wave 
laser and an rf field. The situation is different from the 
CDT because the process of photon emission is dissi- 
pative and, therefore, causes decoherence in the single 
molecule system. Using the renormalization group (RG) 
method for differential equations jjj, we obtain analyti- 
cal expressions for the average waiting time (r) of pho- 
ton emission. Although, the decoherence due to finite 
life time 1 /Y of the excited state of the molecule destroys 
complete localization, in other words, (r) is always finite, 
we find that (r) exhibits resonance behavior as a func- 
tion of control parameter. The most surprising result is 
the existence of a critical value of Y/uj r f above which the 
resonance peaks (i.e. tendency to localize) disappear. It 
turns out, that the transition from coherent localization- 
like behavior to noncoherent and delocalized behavior is 



not smooth. 

Our description of a single molecule system is based on 
the quantum jump approach pj, For a molecule with a 
ground electronic state \g) and an excited electronic state 
|e), the Schrodinger equation is 



^ !*(*)> = 



V g coscj rf t\g) (g\ + -(\g) (e| + |e) (g\) 
(v e cosu) r ft-i~ + Sj |e) (e|j|#(*)). 
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Here S is the laser detuning, V e , g = M ejS ■ E r f, 
H eg ■ E (the Rabi frequency), where /x e g are perma- 
nent dipole moments of the molecule in states |e) and 
\g), and \i eg is the transition dipole, E and E r f stand 
for the amplitudes of the laser and the rf fields respec- 
tively. The Hamiltonian in Eq. Q is the non-Hermitian 
quantum jump Hamiltonian, used in Monte Carlo sim- 
ulations of Ref. 0] to propagate the "wave function" 
|* (t)) = * g (t) \gj + * e (t) |e) 0. In Eq. (JTJ| we used 
the rotating wave approximation for the fast on-resonant 
laser field, but not for the slow off-resonant rf field. Al- 
though, the non-Hermitian Hamiltonian in Eq. Q seems 
simplistic, it describes single molecule experiment very 
well as we demonstrate below. 

The chief quantity of interest is the survival proba- 
bility P (t) = (* (t) \ |* (i)), i.e. the probability of no 
emission event to occur in the time interval between 
and t. The knowledge of Pq (t) enables us to compute 
statistics of quantum jumps from state |e) to state \g), 
i.e. the average waiting time for emission event 0,0] 



(r)= I'™ (\y g (t)\ 2 + \y e (t)\ 2 )dt. 



(2) 



In what follows, we first consider the case of fl <C Y, uj r f, 
and later, the case of Y < -C uj r f. We solve the 
Schrodinger equation Q perturbatively employing the 
renormalizationgroup method @ (cf. recent solution for 
CDT problem 0| with the RG procedure Q). Note that 
the limit of small Rabi frequency (fl uj r f) corresponds 
to the limit of high frequency driving field considered pre- 
viously in the CDT problem We study the effect of 
the non-Hermitian part of the Hamiltonian on localiza- 
tion behavior in the same limit. 
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Introducing scalings t i— > u> r ft, V e , g i— > V eig /uj r f, 
S i ► 5/ujrf, r i — ► r/w. r /,ri i— ► efl/to r f (where 
e <C 1) we pass to the "amplitude'-variables: 
$ g (t) = c g {t)exp{-iV g sint}, * e (t) 
c e (t) exp {— iV e sini — Tt/2}. Thus, for the case of 
zero detuning 5 = 0, we arrive at the set of coupled 
ordinary differential equations: 

fc= — OG 



dBn = _eH^T ^ J 2 (fl 
dt 2 fl ^ r 2 +4/c 2 ' 

A:— — oo 

where subscript "i?" indicates, that we are solving for the 
renormalized values of A and B. Thus, to O (e 3 ) Ar and 
Br satisfy 

A R (t)=^(0)exp — £ fsf^ , 

k A;— — oo / 



|*=-ta*£ ^ ^ +rt/2 J fc (0- (3) 

k— — oo 

where we expanded exponents in series of Bessel func- 
tions of the first kind Jk 0|- The modulation index 
£ = V e — V g is a parameter governing the amplitude of 
the driving rf field. 

To solve Eqs. we proceed with the application of 
the RG approach The naive solution to Eqs. is 
given by a power series in small e: C (t) = (c g ,c e ) — 
C (°) (t) + ec« (t) + O (e 2 ) , provided that (t) are O (1) 
for all t. At large enough times, however, one or several 
o- n ~) (t) (n > 1) may become greater than c' ' (t). Such 
terms, responsible for the break-down of the naive ex- 
pansion are called secular terms. The RG procedure, we 
carry out below, regularizes the naive expansion by iden- 
tifying the secular terms and eliminating them This 
technique enables us to obtain the long time asymptotic 
solution. 

At O (e°) , the solution of Eq. is simply a constant 
vector c(°) = (A, B) T . Substituting this result into the 
next order equations, and integrating, we have 

c^ = -inB ± j fc(c) (i-e--/^)^^, 

k— — oo 

& = E J ^0(e rt/2+M -l)^r 2 - (4) 

k— — oo 

Further application of the naive procedure produces 
terms, which grow with time t, i.e. secular terms 

c « =-^iV J 2 (0 - + NST, 

g 2 ^ fc vs/ r 2 + 4& 2 

fe= — oo 

4 3) = ~B E 4(0^^+ NST, (5) 

where NST stands for nonsecular terms. The following 
RG equations for variables A and B are intended to reg- 
ularize the naive expansion (to order O (e 3 )) 

dAn = _eH^T ^ j2 (0 
dt 2 R ^ r 2 + 4fc 2 ' 

fc= — OO 



B R (t) = B(0)e,J-^El g (6) 

^ k— — oo J 

At time t = we assume that the molecule is in the 
ground state, hence c e (0) = 0, c g (0) = 1. Applying 
these initial conditions and replacing A and B with their 
renormalized values (Eqs. 0) in c^ ^ and in Eq. we 
obtain O (e 2 ) perturbation results for amplitudes c e and 
c 9 . Then, we switch back to the original variables $ e , 9 (i) 
and arrive at 

y g (t) = e -^ sint - ct + 0(e 2 ) , 
e 2 fi 2 00 r 

C = -j- E J fc (0 r 2 + 4fc 2 ' 

k— — oo 

oo 

* e (t) = -i e n e - iV ' 8int - r '/ 2 - ct E Jfc ^ 

A;— — oo 

/ rt/2+ikt _ i) (r - 2ik) 

xi — rdp — - + °^- (7) 

For convenience we set e = 1, keeping in mind that the 
Rabi frequency fl is small. Then, using Eq. 01 we 
calculate the mean waiting time (r) 

fc— — oo 

In Fig. Q]we plot the scaled mean time T(t), given by 
Eq. as a function of the modulation index £ , and com- 
pare this prediction with the results of numerical solution 
of Eq. 0. We observe a good agreement between the 
numerical and analytical results, indicating that the RG 
approach indeed captured the global behavior of ^ e , g - 
We also observe, that for the smallest value of the de- 
cay rate T (T/u> r f in original units), maxima of (r) occur 
close to the zeroes of the Bessel function Jo (£), but with 
increase in T, the maxima of (r) shift and broaden. 

To determine the positions of the extremums of (r) as a 
function of £ (Fig. and to investigate the effect of the 
dissipation on the emission suppression (enhancement) 
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Figure 1: Scaled mean time between emission Y (r) is plotted 
as a function of modulation index £. Solid line represents 
RG prediction of Eq. JHJ, dots correspond to the numerical 
solution of Eq. JJl- Rabi frequency Q, = 0.1. Radiative decay 
rates are T = 0.5; 1.5; 3.0 (from top to bottom). Grid-lines 
indicate zeroes of Jo (£). 
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Figure 3: The shifts (£ - £„) /Cn found from Eq. © (solid 
lines) are compared to those found from numerical solution of 
Eq. (boxes). They are plotted as a function of Y (Y/ui r f 
in original units) for the peaks close to the first (upper curve) 
and the second (lower curve) zeroes of Jo (the first two peaks 
in Fig. [Q. Critical points are marked by stars. 
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Figure 2: Illustration to Eq. Bessel function Jo (solid 

line) plotted as a function of £ together with the right hand 
side of Eq. 10, G (Y, f). Radiative decay rates are (a) r = 1.0, 
(b) T = 2.5. Grid-lines indicate zeroes of Jo (f). 



criterion, we use Eq. JHJ an d set d (r) /d£ = 0. This con- 
dition results in the following transcendental equation, 



MO 



MO 



E 



J k (0 (J fc -i (Q - J fc+ i (Q) 
r 2 + Ak 2 



fc=i 

= G(r,0- (9) 

Solution of Eq. © yields values of £ which minimize and 



maximize (r). These values correspond to enhancement 
and reduction of photon emission rate 1/(t). 

In the limit of Y — > we recover the well-known 
condition, obtained previously in the context of CDT, 

(£n) = (£n are the zeros of the Bessel function J ). 
For small Y, since shifts in ^ with respect to £„ are small, 
we can expand both sides of Eq. ® in Taylor series 
around points £„. Then, the shifts of maxima are given 

by 



r 2 



Jk (£n) (Jk+1 (60 - Jk-1 (^n)) 



r 2 + 4fc 2 



(10) 

For arbitrary values of Y, Eq. ® has to be solved nu- 
merically. In Fig. |2]we illustrate this solution by plotting 
G(r,£) and Jo (£) versus ^ for two values of Y. In Fig. 
Et (r = 1.0) we observe the graphs crossings close to the 
roots of Bessel function Jo. However, as we increase Y to 
2.5, crossings disappear (Fig. 03)- In other words, Eq. 
(jHJl does not have solutions for this value of Y, and first 
two maxima of (r) do not exist. A similar transition also 
happens with the minima of (r). As it follows from Eq. 
© , the minima are found in the vicinity of zeroes of J\ 
for r close to zero, and disappear for sufficiently large Y 
(see Fig. Eb). 

This outcome is an indication of unexpected critical be- 
havior of (t) in r. Maxima (minima) of the waiting time 
(r) do not disappear gradually with the increasing decay 
rate, due to the peaks broadening. Instead, there exist a 
critical value of Y for each extremum, which is equivalent 
to existence of infinite number of critical points. We note, 
that in the vicinity of critical point (£ cr , T cr ), the straight 
forward expansion of Eq. 10 yields £ cr — £ oc (r cr — Y) 13 , 
with (3 = 1/2. Furthermore, as can be concluded from 
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Fig. the neighboring pairs of minima and maxima dis- 
appear at the same values of T cr , that is, maxima and 
minima flow towards each other and annihilate. In Fig. 
El we plot shifts in the positions of maxima with respect 
to £ n given by numerical solution of Eq. @ asa function 
of T. Both curves in Fig. El display excellent agreement 
with the results of numerical solution of the Schrodinger 
equation JQ. 

We remark, that the RG procedure for the case of 
nonzero detuning is performed in a similar fashion and 
results in 
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In this case, maxima of mean waiting time occur close 
to zeroes of Js (£) (S is an integer), if L is not large. 
Otherwise, they shift according to the S — scenario. 

So far we have considered the situation when Rabi ex- 
citation frequency is smaller than L. In a view of single 
molecule experiments carried out in Ref. (ltj j. an oppo- 
site situation Q > L is also important. Assuming, that 
ft, T <C uj r f we now turn to the criterion for the destruc- 
tion of emission for this case, and verify the compatibility 
of the RG results and of our model Hamiltonian with the 
available experimental data. 

Note, that if L is smaller than uj r f, we need 
to modify our derivation. This time, we keep the 
laser detuning S in the ee-matrix element of the 
Hamiltonian in Eq. (QJ, and switch to "ampli- 
tudes": # s (t) = c g (t)exp{-iV g sin UJ r ft}, * e (*) = 
c e (<) exp {—iV e sinidrft — Tt/2 — iSt}. We use different 
scaling to ensure that the rf field frequency is much 
greater than both fl and L (i.e. L i— > eT/uj r f, Q 
efl/wrf). Our system of ordinary differential equations 
now reads 



dc g (t) 



dt 



~2~ 



-ikt — i5t 



dc e (t) 
dt 



E 

Aj— — oc 



ikt-\-iSt 



Jk (0 Ce (t) , 



Jk(0 Cg (t)-—C e (t), 



where all the quantities are in units of uj r f . Unlike the 
previous case (Eqs. (@I EJ), the secular terms appear 
(at integer values of 5) in the naive expansion already at 
the first order. The details of calculations with the RG 
method will be given elsewhere, here we only present the 
final expression for the average waiting time: 



m 2 j k (ty 



r 2 + 2ft 2 j fc (£) z + 4{k-sy 



(12) 



The approximate equality sign in Eq. 11211 indicates, 
that it is only valid when resonances are not overlapping, 
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Figure 4: Comparison of the predictions for (r} _1 (S) given 
by Eq. 11211 (solid lines) with the excited state population cal- 
culated from Bloch equations (Ref. 10]) for various values of 
Rabi frequencies (from bottom to the top): f2 = 0.29; 0.9; 3.2 
in u nits of T. Parameters are taken from experiment of Ref. 
[HJ]: UV//27T = 140 MHz, £ = 1.14, F/2-k = 20 MHz (empty 
triangles). The results are arbitrarily shifted vertically for 
transparency, therefore, units on vertical axis are arbitrary. 
Our predictions are in accord with experiments [TH |. 



which is consistent with our assumption T u r f. We 
see that the criterion for destruction of emission is given 
by Jk (0 = for integer values of the detuning S = k (in 
units of w r /). 

Now we are ready to compare our theory with exper- 
iment of Ref. (lQ| . The experimental results had been 
previously matched with the numerical solution of opti- 
cal Bloch equations We repeat the numerical simu- 
lations carried out in Ref. 0], and compare the results 
with analytical prediction of Eq. ltT2l in Fig. ElllJ- All 
three graphs in Fig. Q]illustrate nearly perfect agreement 
between the theory and experiment, suggesting that the 
observation of critical behavior addressed in this Letter 
may be possible. 

This work was supported by National Science Founda- 
tion award CHE-0344930. 



[1] H.-J. Stockmann. Quantum Chaos: An Introduction 
(Cambridge University Press, Cambridge, 1999), p. 135. 

[2] F. Grossmann, T. Dittrich, P. Jung, and P. Hangi, Phys. 
Rev. Lett 67, 516 (1991); M. Grifoni, P. Hangi, Phys. 
Rep. 304, 229 (1998); R. Bavli and H. Metiu, Phys. Rev. 
Lett 69, 1986 (1992); M. Gitterman, J. Phys. A 34, 10119 
(2001); C. E. Creffield, Phys. Rev. B 67, 165301 (2003). 

[3] D. E. Makarov and H. Metiu, J. Chem. Phys. 115, 5989 
(2001). 

[4] M. Orrit, J. Chem. Phys. 117, 10938 (2002); E. Barkai, 
Y. Jung, and R. Silbey, Ann. Rev. of Phys. Chem. 55, 
457 (2004). 

[5] Ch. Brunei, B. Lounis, Ph. Tamarat, and M. Orrit, Phys. 



5 



Rev. Lett. 83, 2722 (1999). 
[6] L. Y. Chen, N. Goldenfeld, and Y. Oono, Phys. Rev. 
E 54, 376 (1996); T. Kunihiro, Prog. Theor. Phys. 97 

(1997) 179-200; K. Nozaki and Y. Oono, Phys. Rev. E 
63, 046101 (2001). 

[7] M. B. Plenio and P. L. Knight, Rev. Mod. Phys. 70, 101 

(1998) ; J. Dalibard, Y. Castin, and K. Moelmer, Phys. 
Rev. Lett 68, 580 (1992). 



[8] M. Frasca, Phys. Rev. B 68, 165315 (2003). 
[9] M. Abramowitz and I. A. Stegun, Handbook of Mathe- 
matical Functions (Dover, New York, 1965). 
[10] Ch. Brunei, B. Lounis, Ph. Tamarat, and M. Orrit, Phys. 

Rev. Lett. 81, 2679 (1998). 
[11] We also compared numerical solution of Eq. Q with RG 
solution, and found a good agreement. 



